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Abstract 



The focus of the article is on fracture criteria for dynamic crack propagation in elastic 
materials with microstructures. Steady-state propagation of a Mode III semi-infinite crack 
subject to loading applied on the crack surfaces is considered. The micropolar behavior of 
C^ , the material is described by the theory of couple-stress elasticity developed by Koitcr. This 

constitutive model includes the characteristic lengths in bending and torsion, and thus it is 
^_^ , able to account for the underlying microstructures of the material. Both translational and 

y ' micro-rotational inertial terms are included in the balance equations, and the behavior of 

Q \ the solution near to the crack tip is investigated by means of an asymptotic analysis. The 

CJ ' asymptotic fields are used to evaluate the dynamic J-integral for a couple-stress material, and 

the energy release rate is derived by the corresponding conservation law. The propagation 
^vq ' stability is studied according to the energy-based Griffith criterion and the obtained results 

^ , are compared to those derived by the application of the maximum total shear stress criterion. 
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1 Introduction 



In many experimental analyses it has been shown that the mechanical behavior of brittle mate- 
rials such as ceramics, composites, cellular materials, foams, masonry, bones tissues, glassy and 
semicrystalline polymers, is strongly affected by the the microstructure. 

The influence of materials inhomogeneities and defects on the mechanical properties of the 
materials and their interactions with cracks have been extensivel y studied in the framework of 
classical theory o f elasticity by means of hornqgeniz ation theories ( Hashinl . Il959l : lEshelbvl . 119761 : 
Budianskvl . Il965l : iBudianskv and O' Connell 119761 ) . and the effective elastic moduli of bodies 
containing several ei isembles of microstructures have been dete rmined using for instance self- 
consistent methods ( Kachanovl . [19871, 11992 : iHuang et alJ . Il994l ). Furthermore, modern multi- 
scale simu lation approaches have been developed for modelling microstructura l properties of the 
mater ials (JAskes et alJ. I2nn8l. I2nn9l : ISilberschmidtl . I2nn9l : lAndrade et "Zl . 1201 il l . 

In iBigoni and DruganI (|2007l ) it has been shown that classical homogenization results de- 
scribe accurately elastic properties of heterogeneous materials in situations where the bodies are 
subjected to slowly- varying loading and the displacement and stress gradients are small. If high 
gradients are present, standard homogenized materials cannot represent the physical response of 
composite elastic media. Moreover, approaches based on classical elasticity theory cannot always 
predict enough accurately the size effect experimentally observed when the representati ve scale 
of th e deformation field becomes comparable to the length scale of the microstructure (JLakesl . 
19861 ). For example, it has been detected that in presence of stress concentration, such as near 
inclusions and holes, the strength of the material is higher if the grain size is smaller, and that 
the bending a.nd to rsional strength of beams and wires are greater if their cross-section is thinner 



(JFleck et al.l . Il994l ) . Since stress concentration factors derived applying standard homogeniza- 



tion procedures to Cauchy elastic materials depend only on the shape of the inhomogeneities and 
not on its size, they cannot describe accurately these effects. On the other hand, the utilization 
of multi-scale techniques for studying microstructural properties of the materials implies chal- 
lenging numerical computations, and the validation of the results requires a critical comparison 
with analytical solutions and experimental data. 
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Generalized theories of continuum mecha nics, such as m icropolar elasticity (jCosserat and Cosseratl. 

1909 ) , indeterininate couple stress elasticity ( Koiteii . ll964l ) and strain gradient theories ( Fleck and Hutchinson , 



may be considerate as an effort to include rig- 
orously defined characteristic length scales and to study influence of microstructures on the 
material behavior avoiding strong numerical calculations required by multi-scale approaches. 
Indeed, exact analytical formulas for characteristic lengths in couple stress elastic materials 
ha ve been derived via higher order homogenization of heterogeneo us Cauchy elastic inaterial s 
bv iBigoni and DruganI ( 2007! ) and via numerical computations bv JAskes and Aifantid ( 201ll ). 
Analytical solutions derived for couple stress and strain gradient solids can also be used for val- 
idating results of num erical simula tions and analyzing experimental data obtained for materials 
with microstructures ( Laked . Il995l ). 



Indeterminate couple stress elasticity theory developed by JKoiterl (jl964l ) provides two dis- 
tinct characteristic length scales for bending and torsion. Moreover, it includes the effects of the 
microrotational inertia, which can be considered as an additional dynamic length scale. There- 
fore, in order to study crack propagation stability in couple-stress elastic materials, new fracture 



criteri a account i ng for both effects of scale lengths and n iicrorotational ine rtia inus^be formu- 
lated ( Morozovl . 1 1 9841 ). For antiplane crack problems, in iGeorgiadid ( 20031 ) and I Radii ( 20081 ) a 
critical level tc for the maximum total shear stress ahead of the crack tip at which the crack 
starts propagating has been proposed as fracture criterion. This is known as the maximum total 
shear stress criterion, and later in the article we will refer to it as the t^^^ criterion. 

The J-integral for static cr ack problems in couple stress elasticity has been derived by 
Lubarda and Markenscofa ( 2000l ) in the case of an homogeneous material, and bv lPiccolroaz et al 
(J2OI2I ) for an interfacial crack. The energy release rate can be evaluated by means of the con- 
servation of this integral, and further in the paper we will refer to it as ERR. Nevertheless, 
energy-based dynamic fracture criteria for this kind of materials are still unknown in literature. 
For that reason, the principal aim of this paper is to generalize the static energy release rate 
expression to the case of dynamic steady state crack propagation, and to study the effec ts of the 
microstructure on the propagation stability by applying the energy Gr if fith c riterion (JFreund . 



(|20il), 



where the t^ 



19981 ). The results are compared with those obtained in lMishuris et al 
criterion has been adopted. 

The structure of the paper is organized as follows: in Section [2] the problem of a semi- 
infinite Mode III steady state propagating crack in couple stress elastic materials is formulated. 
Both translational and micro-rotational inertial terms are included in the balance equations, 
and a distributed loading app lied on the crack surface s is ass umed . In Sec t ion [3i the dynamic 
conservation laws derived by iFreund and Hutchinson! ( 19851 ) and iFreundl (| 19981 ) for classical 
elasticity are generalized to couple stress materials. General expressions for the dynamic J- 
integral and energy release rate associated to steadily propagating cracks in couple stress elastic 
solids are derived. Explicit forms are obtained for the case of a Mode III crack, and the path- 
independence of the J-integral is demonstrated in Appendix A. 

An asymptotic analysis of the stress and displacement fields near to the crack tip is performed 
in Section [H The contribution of the asymptotic terms to the dynamic energy release rate is 
analyzed in details, the leading term corresponding to finite non-zero energy is individuated and 
an explicit formula for the J-integral evaluated along a circular path surrounding the crack tip is 
derived. The obtained formula involves a constant term depending on the boundary conditions 
of the problem and indicating the amplitude of the leading order term of the asymptotic shear 
stress. This term is evaluated in clos ed form in Secti o n [5] by performing an asymptotic expansion 
of the full-field solution derived in iMishuris et al.l (J2013l ) for the same loading configuration 
applied at crack faces. In this Section, the asymptotics expansion of the full-field solution is 
also used for deriving an alternative equivalent formula for the dynar nic J-int e gral, calculated 
co nsidering the squa r e-shaped path around t h e crack tip introduced bvlFreundl (119981 1 and used 
GeorgiadisI (J2003l ) ; iGourgiotis et al.l (J201ll ) ; lAravas and Giannakopoulod (|2009l ) . The energy 



m 



release rate associated to a steady propagating Mode III crack in couple stress elastic solids is 
compared to the corresponding expression in classical elastic materials. 

In Section [6l the obtained expression for the energy release rate is used fo r studying subsonic 



crack propagation stabi lity. Assuming the energy-based Griffith criterion (jWillisl . Il97ll . 11967 



Obrezanova et al 



2002! ). under the considered loading conditions the steady state propagation 
turns out to be unstable regardless of the va lues the micros t ructu ral parameters. This result 
appears to be in contrast with those detected in lMishuris et al.l (|2013l ) adopting the t™^^ criterion. 



which instead shows a stabilizing effect in presence of relevant microstructures contribution. In 



the authors' opinion, this discrepancy may be due to the fact that the energy release rate depends 
only on the leading term of the asymptotic expansion of the stresses, which dominates very close 
to the crack tip but provides unphysical features such as negative total shear stress ahead of the 
crack tip. Therefore, at a characteristic distance from the crack tip the sole leading order term 
may not describe the correct behavior of stresses and displacements. Differently, the total shear 
stress involved in the t"^^^ criterion is calculated by means of the full-field solution, that takes 
fully into account the microstructural contributions. As a consequence, in order to study the 
crack propagation stability in elastic materials with microstructure, fracture criteria including 
also higher order terms of the asymptotic stresses and involving two or more characteristic 
parameters should be used. Note that in classical elasticity fracture criteria including highe r 
order terms contributions such as T-stress criterion ( Hancock and Dul . ll99ll : ISmith et al.l . l2006l ) 
have also been proposed. 



2 Steady-state cracks in couple stress elastic materials 



In this Section the problem of the steady-state dynamic propagation of a Mode III crack in 
elastic materials with microstructures is formulated by means of the fully dynamical version 
of the couple-stress elastic model, accounting both translational and micro-rotational inertial 
terms into the balance equations. Reference is made to a fixed Cartesian coordinate system 
{0,xi,X2,x-^) centred at the crack tip at the initial time t = 0. Under a ntiplane shear deforma- 
tion, the indeterminate theory of couple stress elasticity ( Koiten . Il964l ) adopted in the present 
study provides the following kinematical compatibility conditions between the out-of-plane dis- 
placement n3, rotation vector (p, strain tensor s and rotation gradient tensor %: 
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Therefore, the rotations are derived from displacement. The rot ation gradient tensor x^ ^-Iso 
known as deformations curvature tensor or torsion- flexure tensor (JKoiterl . Il964l ). is defined in 
the general three-dimensional case as % = Vy? = V x e. The vanishing of the Saint Venant 
tensor (or incompatibility tensor) requires x to be irrotational: 



Vxx = VxVxe = 0. 



(3) 



Using expressions (l2|), it can be immediately verified that relation ([3]) is satisfied. According 
to the indeterminate couple stress theory the non-symmetric Cauchy stress tensor t can be 
decomposed into a symmetric part a and a skew-symmetric part r, namely t = a + t. In 
addition, the couple stress tensor /i is introduced as the work-conjugated quantity of x^ . The 
reduced tractions vector p and the couple stress tractions vector q are defined as 



r 1^ 

p = t n + -Vfinn X n. 



q = M n - /i^ 



,n, 



(4) 



respectively, where n denotes the outward unit normal and //„„ = n • /in. The conditions of 
dynamic equilibrium of forces and moments, taking into consideration rotational inertia, and 
neglecting body forces and body couples, write 

C^13,l + 0-23,2 + Tl3,l + T23,2 = PU3, /^11,1 + Cr2l,2 + 2r23 = Jc^l, /il2,l + Cr22,2 - 2ri3 = Jc^2, (5) 

where p is the mass density and J is the rotational inertia. 

Within the context of small deformations theory, the total strain e and the deformation 
curvature x ^re connected to stress and couple stress through the following isotropic constitutive 
relations 

a = 2Ge + X{tTe)I, p = 2Gf{x' + VX), (6) 

where G is the elastic shear modulus, i and r] the couple stress parameters, with — 1 < 77 < 1. 
Note that for antiplane deformations tre = 0. Both material parameters i and r] depend on 
the microstructure and can be connected to the material characteristic length in bending and in 
torsion, namely 

4 = i/V2, it = ^yiT^. (7) 

Typical experimental values of ih and it for some classes of materials with microstructure can 
be found in iLaked ( 19861 . Il995l ). and analytical expressions for t hese moduli have b e en de rived 



via homogenization of heterogeneous Cauchy elastic materials bv lBigoni and Drugarj ( 20071 ). The 



limit value of r/ = 1 corresponds to vanishing characteristic length in torsion, which is typical of 
polycristalline metals. Moreover, from the definitions ([7]) it follows that It = ib for 77 = 0.5 and 
it = ib = v2 for r] = —1. The constitutive equations of the indeterminate couple stress theory 
do not define the skew-symmetric part r of the total stress tensor t, which instead is determined 
by the equilibrium equations (l5|)2,3. Constitutive equations ([6]) together with the compatibility 
relations ([1]) and ([2]) give stresses and couple stresses in terms of the out of plane displacement 
u-y. 

ai3 = Gu3^i, 0-23 = Gn3,2, (8) 

Wl = -At22 = G£'^{1 + ?7)ti3,12, Ai21 = G£^(W3,22 " ??W3,ll), fJ-U = -^£^(^^3,11 " 1]U3,22)- (9) 

The introduction of ([9]) into ([5])2,3 yields: 

Tl3 = ^An3,i + -ii3,i, r23 = ^Au3,2 + -ii3,2, (10) 



where A denotes the Laplace operator. By means of (^ and (jlOp . the equation of motion ([5])i 
becomes 

f n J 

AU3 - -AAU3 = ^ils - ^^^3- (11) 

We assume that the crack propagates with a constant velocity v along the xi-axis. In this 
case it is convenient to introduce a moving framework x = xi — vt,y = X2, z = X3, and the out 
of the plane displacement can be assumed in the form: 

U3{xi,X2,t) = w{x,y). (12) 



It follows that the time derivative of the displacement w can be written in terms of the 
derivative with respect to x, namely w = —vw^x and thus iv = v^w^xx- Therefore the equation 



of motion (jlip under steady-state conditions becomes: 



Au; AAw = m {w^xx — h^f. ^w^xx) 



(13) 



where m = v/cg is the normalized crack velocity, Cg = yG/p is the shear wave speed for 
classical elastic m aterials, the ch a racter istic length h is defined as h = v/cj) with cj) = yMZ^/J, 
and ho = h/i (see iMishuris et al.1 (J2013l )). 

According to (H]), the non- vanishing components of the reduced traction and couple stress 
traction vectors along the crack line y = can be written as 



P3 = ^23 + 2/"22,x, qi = /^21, 



(14) 



respectively. By using ^2, (El) 1,21 (fTO]) '> and (fH|) . the skew-symmetry of the Mode III crack 
problem requires ahead of the crack tip: 



w = 0, w,yy — rjw^xx = 0, for x > 0, y = 0. 



(15) 



On the crack surface, vanishing of the reduced traction and couple stress traction yield to 
the following boundary conditions for the function w: 



'^,y~^[{'^ + '^~'^^'^^l)'^,xx+w,yy\y = -YiT{x) w^yy - r]w^xx = 0, for x<0,y = 0, 

(16) 
where t(x) is the loading applied on the crack faces, which is assumed to have the following 
form: 

r(x) = ^e^/^, x<0. (17) 

Note that although we discuss here only a specific loading condition, the main conclusions 
reported in this paper have been confirmed for other types of loading. 



3 Dynamic energy release rate 



In thi s Se ction, the dynamic cons e rvatio n laws obtained for linear elastic media by iFreund 



(j 19981 ) and iFreund and HutchinsonI (j 19851 ) are generalized to couple stress elastic materials. An 
explicit integral expression for the dynamic energy release rate associated to steady state cracks 
propagation in elastic solids with presence of couple stress is derived. 

The energy release rate for a dynamic crack has been defined by iFreundl ( 19981 ) as the 
following limit: 

-F{ry 



£ = lim 



(18) 



where F is the total energy flux through the contour T surrounding the crack tip. For couple 
stress elastic materials, the energy flux for a dynamic crack propagating along xi-axis is given 
by: 

{W + T)vni + t^n • -^ + fJn ■ ^ 



F(r) 



ds, 



dt ' ' dt 

where n is an outward unit normal onT, W denotes the strain-energy density 



W 



and T is the kinetic energy density 
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(19) 



(20) 



(21) 



Since a Mode III steady-state crack propagating at constant velocity v along the xi-axis is 
considered, the expression (|12p for the out-of-plane displacement is used and then in the moving 
framework {x, y, z) the energy flux ()19p assumes the special form 



F{T) =v j [{W + T)n.^ - t^n • e,w,^ - fi^n ■ <p.^ 



ds, 



(22) 



then the generalized J-integral for an antiplane dynamic steady-state crack in couple-stress 
elastic materials can be defined: 



J = ffil = [ [{W + T)n^ - t^n • e,w,^ - /i^n • (^ ^.] ds 

= I [{W + r)n - p • e^Vw - (V<^)^q] • e^ds 



(23) 
(24) 



The expressions ()23p and has been proved to be path independent, the details of the demon- 
stration are reported in Appendix A. Moreover, the equivalence of the two alternative forms of 
the J-integral (j23p and ()24p , the first written in function of the tractions and the second of the 
reduced tractions, is de monstrated in Appendix B. The (1231) i s the generalization of the static 



Lubarda and Markenscofi 



expres sions derived by lAtkinson and LeppingtonI (jl974l . 119771 ) , and 

(2000) to the antiplane dynamic steady state case. The dynamic energy release rate is then 
defined by the limit: 



£ = lim / \{W + T)nx — t'^n • GzW x — At'^n • (p J ds. 

Using definitions ([T|) and ([2]), the strain energy density (j20p becomes 

G£2 



W 



Cx 2 2 

"^l^,x + ^,j/; "I Y' LV"^,a;a; T i^,yy) -r ^\^ -r 'l)K"J,xy ^ "^,xx^,yy 
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whereas for steady state propagation the kinetic energy density is given by 



T= — 
2 



PW% + j{W%y + W%) 



(25) 



(26) 



(27) 



A polar coordinates system (r, 9) centered at the crack tip is assumed, and a circular contour 
of radius r around the crack tip with n = (cos ^, sin 0,0) is considered. Then, substituting 
expressions ([26]) and ([27|) the J-integral (|23|) becomes 



J 
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pw^^ cos 6 + —{w^xy + w^^xx) COS 6 - ^w^^ iw,xxx COS 6 + w^y^x sm9) 



J_ 
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2w^xW^y sin 0] + 



(28) 



ce 



(Aw)' 



cos 9 — Aw {w^xx cos 9 + w^xy sin 0) + tf^x {Aw^x cos + Aif^y sin 9) 



rd9. 



According to the definition ()25p . the energy release rate can be evaluated as the limit for 
r — )• of the integral ([28]). In order to evaluate explicitly the J-integral ([28p and to investigate 
the variation of the energy release rate ([25p in function of the crack propagation velocity and 
its implications on propagation stability, the behavior of the out-of-plane displacement w near 
to the crack tip is studied by means of an asymptotic analysis in the next Section. 



4 Asymptotic crack tip fields 

The following standard asymptotic expression for out-of-plane displacement w in separate vari- 
ables form is considered: 

w{r,9) = r'Fs{9), r ^ 0, (29) 

We are interested in finding the terms of the asymptotic solution ([29p corresponding to finite 
and non-zero contributions to the J-integral ([28p in the limit r — ?• 0. Since the displacement 
w should be bounded and symmetrical, it follows immediately that s > 0, and then no zero- 
order terms are present in the asymptotic expansion. Substituting the expression (|29p in the 
generalized J-integral formula ([28p and using the following derivative rules which hold for an 
arbitrary function f{x,y) = f{r,9): 



f,x = f,r cos 9 - /_, 



sm( 



f,y = f,r sin 9 + f^e 



cost 



(30) 



we get: 
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2s-l f-K 
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^2s — 3 fTT 



[pv'^{sFs cos 9-F'^ sin 9f cos9 + G (1 + s)F,F,' sin 9 - {sF^ + f','^) cos ^1 1 d^ + 



4 



+ ^— — / {of {s'^Fs + f',' )\{Fs-s{2- s)f'^) cos 9 + 2F', sin 9 



+ 



+ — -(2 - s){sFs cos 9- f'^ sin 61) s{s - (2 - s) cos29)Fs + 2(1 - s)f', sin 26* + 2FJ sin^ i 

where the superscript denotes the total derivative respect to the vari able 9 . 

Observing expression ([^T]) . in agreement with the results reported in lRadil (|2008l ) for the static 
case, we deduce that the finiteness of the energy release rate towards the crack tip requires that 
s > 3/2 for any non integer number. 



(31) 
d9, 



Including higher order terms in the asymptotic expansion in the form ^- r'^'-Fs-{9) and using 



this expression in ()28p . terms of order r^i+^j-^ and r^^'^^^~'^ where i ^ j are detected. These 
terms involve both Fg^ and Fs ■ and can have a non vanishing impact to the value of the energy 
release rate. Therefore we need to consider several asymptotic terms in the form of (j29p and 
analyze the possible correlation between them in the nonlinear functional ()28p . According with 
this discussion and in order to find the terms corresponding to finite and non-zero contributions 
to the energy release rate, we assume 1 < s < 3. 

It can be easily demonstrated that if |sj — Sj| < 2, V i 7^ j, only the leading order term 
of the governing equation (|13p can be considered, while if more terms are required with expo- 
nents differing by 2 or more than 2, the full equation ()13p must be considered in the analysis 
( Piccolroaz et al.l . I2OI2I ). As a consequence, assuming 1 < s < 3, we can then keep only the 
leading term of the evolution equation (fT3|) 



A{Aw 



X^w^ 







where 



A^ 



2m^hl 



2GP' 



(32) 



Introducing the expression (|29p in (|32p . the general asymptotic solution of the equation of 
motion has been obtained, the derivation is illustrated in details in the Section containing the 
supplementary material. Referring to this general solution, since we are assuming values in the 
range 1 < s < 3, the terms corresponding to s = 1,3/2,2,5/2 needs to be considered and the 
asymptotic expression for the out of plane displacement w turns out to be: 



w{r,9) = Birsme + B2r^^'^ 



sin h-{l- A2 sin^ ^)3/4^±l^ sin ^$(0) 



sm ■ 



:i-A^sin2 0)V^-l±^sin|$(0) 
1 + ry — A^ z 



+ (33) 



where 



^{6) 



arcsm 



^/r^A2 



smi 



Vl - A2 sin2 9 



(34) 



and the constants Bi,B2, B^ and B4 define the amplitude of each asymptotic term the sum and 
thus must to be evaluated according to the boundary conditions (|15p and (J16p . 

It has been verified by symbolic calculations tha t the terms o f the (1331) corresponding to s = 1 
and s = 2, similarly to what has been detected bv iRadil (|2008l ) and IPiccolroaz et al.l (|2012l ) for 
a stationary crack case do not contribute to the J-integral and to energy release rate. Although 
these terms are relevant for evaluating displacement and total shear stress at the crack tip, the 
only finite and non-vanishing contribution to the generalized J-integral (|28p is associated to the 
order s = 3/2 of the asymptotic expression (j33p . and then the energy release rate is given by: 
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(35) 
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6*) 3(1 - 3 cos 26l)i:^3/2 - 4iJ3/2 sin 29 + SHy^ sin^ 9 | 
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where 



Hy,i9) = sin ^^ - (1 - X' sin2 0)3/4 _^±^ sin ^a>(0) = (36) 



sin-0-^ f^l-A2sin2 + 2cose') V Vl - A2 sin^ - cos( 

2 y2fi + ?7-A2) V^ ; V V 



The proposed procedure, based on the assumption 1 < s < 3 and then on the analysis of 
the leading order term of the governing equation (j32p . can in principle be performed considering 
a different range of values of s, as 1 < s < 3, corresponding to terms s = 3/2,2,5/2,3, but 
observing the expression (j3ip and remembering that the other contributions to the J-integral 
are of orders r^^~^^^~^ and r**"''*^"^ ^nd that s = 1/2 is excluded because finite energy is required, 
it is easy to deduce that all terms associated to s > 5/2 do not contribute to the energy release 
rate. Therefore we can conclude that s = 3/2 provides effectively the only term contributing to 
the J-integral and that the choice of considering the leading term of the governing equation (|13p 
in the asymptotic analysis is correct for our purpose. 

In order to evaluate the energy release rate ()35p . the constant B2 must be explicitly de- 
termined. In the next Secti on B2 is calculated s tarting from the asymptotic expansion of the 
full- field solution derived in lMishuris et al.l ( 20131 ) for the same loading conditions (fT5]) -(fT7 |) by 



means of Wiener-Hopf technique. 

5 Explicit evaluation of the energy release rate 

In order to derive an explicit expression for the constant B2, we perform the asymp totic analysis 
of the Fourier transform of the full-field solution derived in iMishuris et al.l ( 20131 ) for the same 



loading conditions ()15p -(fT7|) by means of Wiener-Hopf technique. In the limit |s| — t- 00, the 
Fourier transforms of stress, couple stress field and displacements assume the following behavior: 

tUs, 0*) = ^^°^<' + " - ^f "'> (sO'f + O {(si)-"'') . Im,,>0. (37) 

T(/io,m,?7) ^ V ^ / 



2iFTo^HU/l- 2/^2^2-7? (1 + 7?) 

71+ (5,0+) = ^ ^ , {si)-'^' + 0{is£)^'), lms>0, (38) 

T(/io, m,7])(l + ^l-2^m? 



2FT0G , 5/2 , ^ /, ,,-7/2 



^"(^,0+) = -— — ^^ r(5^)r^^ + (5^):'^^ , lm.<0, (39) 

Gl [no, m,rjj \ / 

where F is a constant determined starting from the the boundary condition ()15p i and applying 
the Liouville theorem and: 



K{i£/L) 1 - 7?2 - 2hlm^ + 2v/l - 2/^2^2(1 + r?2 - 2hlm^) 

{ii/Lyr 1 + Vl - 2/i2m2 
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the explicit expression for the factorization function k^{ii/L) is given bv lMishuris et alJ ( 20131 ). 
Substituting the <^(i) into ([37]) and ([39]) we obtain: 



2iFTok+iii/L) Ul - 2/igm2 - r? (1 + ,?) 
71+ (s,0+) = ^ ^ ' ^ is£)-'^' + O (s£);i , Im. > 0,(42) 



uI-(.,0+) = ?^^^§±(^^/^L(.£)-^/V0 f(.£)~^/^) , Im. < 0. (43) 

Giii/L)X^Tiho,m,7j) ^ ' 



Further, we consider the fohowing transformation formula ( Rood . Il969l ): 



^K /* ^K+ip(^ + l)s-''-\ with K / -1, -2, -3 . . . , (44) 

where F is the gamma function and the symbol -H- indicates that the quantities on the two 
sides of the ()44p are connected by means of unilateral Fourier transform. Applying the ()44p to 
expressions (j^Tj) and ([15|) . we get: 

*,3(..0-) = -£I°^|^±(!^<l;;i^^.-'^ .>0. (45) 

2/7r T{ho,m,ri) 



2FToVLk+{zi/L) Vl - 2/igm2 - r? (1 + t?) 

//22(x,0+) = ^^ ^ , x-^/^, x>0, (46) 

V^T(/io, m, r/) 1 + ^1 - 2/i^m2 



On the crack surface, for = tt and y = O"*", the term of order 3/2 of the asymptotic 
expression in polar coordinates ()33p becomes 



^^'^'°^^ = ^^rT7^5^^"'^^'^'' '^^°' ^^^^ 



where the definition r(y = 0+) = v j;^ = \x\ has been used. Equating expressions (j^7|) and (08]), 
we get: 

_ 4FTo^/Lfc+(i£/L) / l + r?-2m2/.g \ 
^ 3^GT{ho,m,r])i^ \ m?hl )' ^ ' 



The explicit expression (j49p can the be used into the (j35p for studying the variation of the 
energy release rate in function of the crack tip speed and of the microstructural parameters /iq 
and rj. 

In order to check the validity of the obtained results, an alternative expression for the energy 
release rate is derived considering the rectangular- shaped contou r F with vanishing height along 
the y-direction and with e — )• +0 reported in FigJT]( Freund . ll998l ). The Cartesian components of 
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the outward unit vector normal to T are n = {nx,ny,0). Considering the moving framework in 
Fig{T]with the origin at the crack tip, for the steady-state antiplane crack problem the generalized 
J-integral (f23|) becomes: 



J = [{W + T)nx - (tisnx + t23"-j/)«^,x- - {fJ-iinx + Ai2i%)v?i,z - (/^i2"-z + /^22%)9'2,x-] ds = 

= / [{W + T) -tisW^x- {flinpi^r,+ Hi2(p2,x)]dy - / [t23W,x + {fJ'2l'Pl,x + ^J'22V^2,x)]dx,{50) 

where ti3 and ^23 are components of the total non-symmetric Cauchy stress tensor, including 
both symmetric and skew-symmetric part. 
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y '^ 




(e,+0) 
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(-B.-0) 



(e,-0) 



Figure 1: Rectangular-shaped contour around the crack-tip 

In order to evaluate the energy release rate, we allow the height of the rectangular path 
reported in Fig{T]to vanish. In this limit, the first integral of the ()50p is zero and the following 
expression for the dynamic energy release rate is derived: 



£ = -2 lim <^ / [t23W^x + {fJ'2i'Pi,x + fJ'22'P2,x)] dx \ . 



(51) 



It is important to note that anti-symmetry conditions (jlSp together with boundary conditions 
PU|) provide that the reduced traction qi = ^121 is zero along the whole crack line y = 0, where 
n = (0, ±1, 0). Consequently, the energy release rate ([5T]l becomes: 



£ 



-2 lim 



-2 lim 

e->+0 



-2 lim 



[t23(2;,0+)i(;,a;(x,0+) -F/U2i(a:;,0+)(^i,x(3;,0+) + /i22(a:,0+)v32,x(2;,0+)] dx 



[t23{x, 0+)ti;,^.(x, 0+) + ^jL22{x, 0+)992,:,.(2;, 0+)] dx 

t2z{x, 0+)tf,a;(x, 0+) - -fJ.22{x, 0^)w^xxix, 0+) dx 



(52) 
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For evaluating the integral ()52p . solely asymptotic expressions for the traction ahead of 
the crack tip ^23, the couple stress field 1122, and the displacement w are required. Then, by 
substituting equations ([15]) . (P6]l and ()17|l in the general formula ([52]) . we finally derive: 



S 



-2 lim 



2F^T^Lkl{ii/L) 
' e^+o [ TrGi^T^{ho,m,T]) 



(1 + V 



2hlrr?) 



+e 



1/2 -3/2 , 



1 - him? + 1 



X x,dx 



(53) 



where x_ ,x_ and x, jX, are distributions of the bisection type ( Arfken and Weber 



2OO5I ). For any real a with the exception of a = 1, 2, 3, 
is defined as follows: 



X, 



|x|", for X > 0, 
0, for X < 0. 



. , this particular type of distributions 

0, for X > 0, 

Ixl", for X < 0. 



Moreover, the product of dist ributioiis insid e the integral in ([53|) is evaluated through the ap- 
plication of Fisher's theorem ( Fischeii . Il97ll ). that leads to the operational relation: 



(x_)"(x+) 



-1-0 



it6{x) 
2sin(7ra) 



with a 7^ —1, —2, —3 . . . , 



(54) 



where 5{x) is the Dirac delta distribution. Then, by using the relation ()54p int o (1531) and consider- 
ing th e fundamental property of the Dirac delta distribution J_ ^ 6{x)dx = 1 ( Arfken and Weber 
200i), we finally get: 



£ 



2F'^T^Lkl{i£/L) / 2(1 + r/ - /igm^)^! - 2/i^7ti2 + (1 - 2hln? 
G£2T2(/io, m,r/) I 



2/i2m2 + 1 



^\ _ 2F'T^Lki{ie/L) 
~) Ge^T{ho,m,i]) ' 
(55) 
Using this alternative procedure, we have derived the explicit expression ([5^ for the energy 



VT 



release rate corresponding to a Mode III steady state propagating crack in a couple stress elastic 
material. Equation (|55p can be compared with the energy release rate associate to an antiplane 
steady state crack in a classical elastic material, derived assuming the same loading configuration 



£• 



d 



rp2 



1 



GL^T 



m^ 



the ratio between the two expressions (|55p and ()56p is given by: 

£ 2F^L'^kl{U/L) 



£^i PT{ho,m,ri) 



■v^ 



m^ 



(56) 



(57) 



The analytical expression (j55p has been proved to be equivalent to (|35p by means of several 
numerical examples. Indeed, both expressions ()35p and (|55p provide the same results, which are 
reported in the next Section. 
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6 Results and discussion 



In this Se ction, the var iation of the energy release rate is analyzed applying the classical Griffith 
criterion (jWillisl . 119671 ) in order to study the propagati on stability. The resu lts are compared 
with those detected using the t'^'^^ criterion, adopted bv lMishuris et al.l (J2013l ). 

The normalized variation of the energy release rate versus the crack tip speed m is reported 
in Fig. [2] for the same value of the ratio L/i = 10, three different values of rj = {0, 0.9, —0.9} 
and four different values of the rotational inertia /iq = {0.01, 0.6, 0.707, 0.8}. The range of the 
normalized crack tip velocity has been chosen in such a way t hat the propagation is subsonic and 
the conditions of validity of the full-field solution obtained in lMishuris et al.l ( 20131 ) and used for 
evaluating the J-integral constant are satisfied. A similar behavior is observed for all different 
set of parameters: the energy release rate is initially constant for small values of the crack tip 
speed m < 0.3, then increases monotonically until its limiting value. For small values of /iq and 
r/, the limit value of the energy release rate corresponds to ttt. = 1, while as the microstructural 
parameters increases the limit value of the crack tip speed becomes smaller than the shear waves 
speed Cs, and thus m < 1. 
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Figure 2: Variation of the energy release rate as a function of tiie normalized crack tip speed 
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Figure 3: Variation of the ratio £ /£'^^ as a function of the normalized crack tip speed m. 
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On the basis of the Griffith criterion ( Freundl . ll998l ). crack initiation requires that the energ; 



release rate achieves a critical value Ec = 27, where 7 is the energy needed to form a unit of 
new material surface, and is a constant depending only by the properties of the medium. In our 
case, once this critical value for the crack initiation is achieved, the energy release rate always 
increases in function of the velocity, this means that if the applied loading provides the energy 
necessary for starting the fracture proce ss, the crack has enough energy to ac celerate rapidly 
up to the limiting values of the speed ( Willid . Il97ll : lObrezanova et al.1 . |2002| ). This implies 



that, if we analyze the results shown in Fig. [2] by means of Griffith criterion, the steady state 
propagation of a Mode III crack in couple-stress elastic material can be considered unstable for 
any value of the rotational inertia /ig and of rj. 

The ratio between the energy release rate in couple stress materials and the energy release 
rate in classical elastic materials ()57p is plotted in Fig. [3] as a function of the normalized crack 
tip speed m. For small values of the rotational inertia, this ratio decreases as m tends to 1. This 
is due to the fact that for m = 1, the energy release rate corresponding to a classical elastic 
material ()56p diverges, while in presence of couple stress it has a finite limiting value. Differently, 
for large values of /iq, the ratio increases monotonically until a limiting value corresponding to 
m < 1 and thus to a crack tip speed smaller than c^. Therefore, in order to propagate the 
crack at constant velocity, more energy than in a classical elastic material must be provided by 
the loading if the contribution of the rotational inertia is relevant. Observing Fig. [3l we can 
also note that in the static limit, namely as m tends to zero, the ratio £ JS'^^ is close to one for 
r] = —1, this behavior is in agreement with the results reported bv iRadil ( 20081 ) . which illustrate 
that f or 7? = — 1 the s o lution approaches the cla ssical elastic case. 



In iMishuris et al.l ( 20131 ) the t™^^ criterion ( Georgiadisl . |2003|) is applied to the same crack 



problem. This alternative criterion states that the maximum total shear stress ^23^^^ possesses 
a critical level tq at which the crack starts propagating. The behavior of ^23^^ ^^ ^ function of 
the crack tip speed has been studied extensively for several sets of microstructural parameters. 
For each values of the ratio rj, it has been individuated a critical value of the rotational inertia 
hoc, such that for ho > /iqc the maximum shear stress increases very rapidly in function of m 
and becomes unbounded at a crack-tip speed lower than the shear wave speed Cg. In these cases 
the crack propagation turns out to be unstable and a limit speed of propagation Vc < Cg is 
individuated. For ho < /iqc, ^23^^ decreases instead as the crack speed increases and tends to 
zero as m approaches its limiting value, suggesting the occurrence of stable crack propagation at 
velocities sufficiently lower than the shear wave speed. The critical value of the rotational inertia 
is reported as a function of rj in FigH] (blue line). The limit value 77 = —1 corresponds to hoc = 0, 
then in this case the propagation turns out to be unstable for any value of the rotational inertia, 
while as r] increases hoc grows, and the range of the rotational inertia associated to stable cracks 
propagation becomes larger. Therefore, the present analysis shows that if the maximum total 
shear stress criterion is adopted, a stabilizing effect of the crack propagation is provided as the 
characteristic ratio rj increases and then the contribution of the microstructures effect becomes 
relevant. Differently, for negative values of rj near to the limit case r/ = — 1, crack propagation 
is detected to be unstable as the rotational inertia becomes not negligible. 

The stabilizing effect is not detected applying the Griffith criterion, according to which the 
propagation is unstable regardless of the values of rj and Kq. As a consequence, assuming this 
criterion, the critical value of the rotational inertia is zero for all values of rj (dashed line in 
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Fig. HI), and the fracture is unstable for each set of microstructural parameters. This is due to 
the fact that the energy release rate is evaluated using the term of order r^'"^ of the asymptotic 
displacement field, corresponding to the singular shear stress t erm of order r~^/ ^ . As i t has been 



illust rated by many studie s both in class i cal elastic materi als (jDu and Hancockl . 1 19911 : iDu et al. 



199l|) and in couple stress ( Radi and Geil2004l : lRadil . l2008l ). this singular contribution dominates 
very near to the crack tip, but it displays an unphysical negative shear stress ahead of the crack 
tip. Moreover, it is not sufficient to describe accurately the physical behavior of the stresses 
at a characteristic length from the crack tip where the higher order terms of the expansions 
become important. In particular, the second term, which is associated to the linear term of the 
displacement and is commonly known as T-stress, can infl uence significantly the processes o f 
crack initiation and propagation in many physical situations ( Tvergaard and Hutchinsonlll994l ). 
In th is cases, the c ritical stress intensity factor criterion and thus the connected Grifflth crite- 
rion ( Freundl . Il998l ) are insufficient for describing accurately the crack extension and stability. 
Alternative two-parameters fracture criteria requiring the achievement of a critical value for 
both the stress intensity factor and the T-st ress has been proposed for classical elastic materials 
(JHancock and Dul . Il99ll : ISmith et al.l . l2006l ) and can be extended to couple stress. 



0.8 
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Figure 4: Critical value of the rotational inertia /loc as a function of the ratio 77. 



In couple stress materials, as the distance from the crack tip increases the discrepancy be- 
tween the stresses physical behavior and the sing ular l e ading term of the asymptotic is even 
more relevant respect to the classical elastic case ( Radii . 12008 ). Therefore, the contribution of 
the higher order term in the asymptotic shear stre ss, which does not affect the energy release 
rate, is relevant. Differently, in lMishuris et al.l ( 20131 ) the total shear stress is evaluated by means 
of the full-field solution, that describes correctly the behavior of displacement and stresses on 
the whole crack line and takes fully into account the action of the microstructures. 
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Fig. [5] shows that for a fixed value of the crack tip speed, the energy release rate slightly 
increases its magnitude for large values of Kq. As a consequence, for large values of the rotational 
inertia, a greater flux of energy at the crack tip is detected, and thus crack initiation and propa- 
gation are favored. Moreover, in Fig. [5]we can also observe that ERR and the ratio E jEc decrease 
as T] increase. In particular, the ratio E jEc is almost one for 77 = —1, and then it decreases. As 
it is reported in Fig. [6l this behavior is observed also in the stati c case, corresponding to the 
limit TTi = 0. This confirms the results illustrated in iRadil ( 20081 ) and can be explained with 
the fact that for large values of ry, associated with large values of the characteristic length in 
torsion, the toughness of the material increases, and thus less energy is provided at the crack 
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tip and a shielding effect due to microstructures is detected. It is important to note that in the 
static case the energy release rate is independent by the value of /iq , because the rotational iner- 
tia appears only in the dynamical terms of the evolution equation ()13p and of the J-integral ()28p . 



7 Conclusions 



A general expression for the J-integral associated to dynamic steady-state cracks subjected 
to antiplane loading in couple-stress elastic materials has been derived. The effects of both 
finite characteristics lengths as well as rotational inertia are included in the analysis. The 
generalized J-integral has been demonstrated to be path independent for steadily propagating 
cracks. Asymptotic displacements and stress fields have been derived and used for evaluating 
the dynamic energy release rate. Surprisingly, the dependence of the energy release rate by the 
microstructural parameters looks less than the authors expectations. However, for large values 
of the characteristic length in torsion the ERR decreases, and thus less energy is provided at the 
crack tip and a shielding effect due to microstructures is detected. 

The s tability of the propagation has been stud ied by applying the energy-based Griffith 
criterion ( Willid . 119671 . ll97ll : IObrezanova et al.l . l2002l ). The propagation turns out to be unstable 
for each value of the characteristic length i n torsion and o f the r otational inertia. This result 
appears to contrast with those detected in iMishuris et al.l ( 20131 ) . where the maximum total 
shear stress criterion has been adopted, and a stabilizing effect in correspondence of relevant 
microstructural contributions is shown. In the authors' opinion, the discrepancy is due to the fact 
that the energy release rate depends only on a single term of the asymptotic displacement field, 
corresponding to the singular leading order term of the stresses. The singular stress dominates 
in proximity of the crack tip, whereas as the distance from the crack tip increases the effects of 
higher order asymptotic terms which do not contribute to the J-integral become relevant. It is 
important to note that analogous re sults has been detected by compar ing the ERR and the T- 
stress criterion in classical elasticity ( Du et al.l . ll99ll : ISmith et al.l . l2006l ). Moreover, the presence 
of the microstructures affects stresses and disp lacernent up to a distance from the crack tip of the 
order of 5-10 times the characteristic length i ( Radi 120081 ) . and thus cannot be described by the 
sole leading order term. Differently, the total shear stress is calculated by means of the full-field 
solution, that takes fully into account the microstructural contributions. Therefore, in order 
to study cracks propagation stability in elastic media with microstructures, fracture criteria 
considering also higher order terms must be developed and validated by means of experimental 
analysis. The t™^^ criterion may represent one of the possible alternative methods that could 
give relevant results concerning microstructural effects on crack propagation stability. 
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Supplementary Material 

In this Section a general asymptotic solution for the evolution equation ([32]) is derived. Equation 
(|32p can be split in the following two PDEs 

Aw = 0, Au; - \^w^^^ = 0. (58) 

Substituting the asymptotic expression for the displacement (|29p into ()58p i . and using the deriva- 
tive rules (j30p . the following ODE for the unknown function Fs{9) is derived: 

f'^{9) + s'Fs{9) = ^. (59) 

This equation admits the solutions providing: 

w\r, 6) = r'[Ci cos{s9) + C2 sm{s9)], (60) 

for any s > 0. 

Equation ([58]) 2 can be reduced in the form ()58P i by considering the linear coordinate trans- 
formation 

X = x, Y = yVl- A2. (61) 

Let us denote the transformed coordinates as X = Rcos $ and Y = i?sin$, where: 



R(r, 9) = \/X'^ + y2 = ^^2 + (1 - A2)y2 = j.\J\ - A2 sin2 9, (62) 

(Y\ ( Vl-A2sin0 \ , , 

<P(y) = arctan -— = arcsin — . (63) 

\X) VVl-A2sin2 0y/ 

Then, the general solution of the equation (|58p 9 is given by 

u;^^(r,^) = i?^[C3Cos(s$(^)) + C74sin(s$(0))], (64) 

for any s > 0. 

Finally, the general solution of the governing equation ()32p is obtained by the sum of w and 
w^^ ^ given respectively by (|6Up and (|64p . By using expression (|63p . we get: 

w(r, 9) = r'[Ci cos{s9) + C2 sin(s0) + (1 - A2 sin2 9y/^{C3 cos{s^{9)) + C4 sin(s$(0)))]. (65) 

Since the asymptotic expansion of the loading function (|17p contains only integer powers of 
r, for non-integer values of the exponent s, which in our case means for s 7^ 1,2, the boundary 
conditions ()15p at 9 = and (|16p at 6 = ir yield the following linear and homogeneous system 
for the u nknown constants C\,C2,C3 an d C4, equivalent to that derived for the traction- free 



problem (JRadi and Geil . l2004l : iRadil . l2008l ^: 



Ci + C3 = 0, 

sis - 1)[(1 + r?)Ci + (1 + r? - A2)C3] = 0, 

sis - 1) { [(1 + r])Ci + (1 + ?? - A2)C3] cos(s^) + [(1 + 7])C2 + (1 + r? - A2)C4] sin(s7r)} = 0, 

sis - l)(s - 2) { [(1 + r/ - A2)C2 + VT^J?il + r])Ci] cos(s7r) 

-[(1 + r/ - A2)C7i + VT^X^il + r/)C73] sin(s7r)} = 0. 

(66) 
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Then from ()66p it necessarily follows that Ci = C3 = and 

[(1 + T])C2 + {l + v- A2)C4] sin(s7r) = 0, 



[(1 + 7? - \^)C2 + \/r^A2(l + 7?)C4] cos(s7r) = 0. 



(67) 



Eqs. (|67p admit a non trivial solution for sin27rs = 0, namely for s = njl where n € N. In 
particular, for n odd one obtains from (j67P i : 



C4 



1 + r? 



1 + 7? - A2 



C2, 



and consequently from (|65p 
u;(r, ^) = C2r"/2 

For n even, from (|67p 9 it follows: 



n 
sm I —( 
2 



C4 



1+r/ 



1 + r? - A2 



(1 - A^ sin 



2flW4. 



n _ 



sin -cl>(0) 



1 + 77 - A^ 



and consequently from (j65 



1 + r/ - A^ 



C2, 



sm ( — „ , , 

2 ; vr^A^(i + ^) 



(1-A2sin2 0)"/Sin -$(0) 



n . 



(68) 



(69) 



(70) 



(71) 



Since we are considering values in the range 1 < s < 3, the terms corresponding to s = 
1,3/2,2,5/2 need to be included in the asymptotic solution. The terms of the order s = 3/2 
and s = 5/2 possess the form ([69]) . while the terms s = 1,2 correspond to degenerate cases of 
the equation ([32]) . and need to be treated separately. 

For s = 1, the general solution of the equation (j32p may be found by solving the problem 



^w = r~i(l - A^ sin^ ey^l^{C^ cos ^{6) - Q sin $(^)), 



namely 



Auj = r^i (1 - A^ sin^ 6l)~^(C3 cos 6 - C4\/l-A2sin ( 



(72) 



(73) 



By using the method of variation of the parameters, one can find the following solution of 
equation ([73]) in the separate variable form (p9]) 



w[r, 



r 
2A2 



Ci + C3 log(l - A^ sin^ e) + 2Ci{^ - 0\/l - A^) 



cos 9+ 



+ 



C2 - 2C73($Vl - A2 -9) + C4\/l-A2log(l - A^ sin^ 



sm( 



(74) 



The boundary conditions ()15p and ()16p necessarily imply that Ci = C3 = C4 = 0, so that 
the corresponding solution reduces to 



w[r. 



Car sin 6*, 



(75) 
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where C2 depends by the amphtude of the loading (|17p and also by the constant associated to 
the term s = 3, that is neglected in our analysis. 

For s = 2, the general solution of the equation ()32p may be found by solving the problem 



Au' = C3 + C4$(0), (76) 

namely 

r. r. I \/r^^A2sin0 \ 

Ai(; = C3 + Ca arcsm = , (77) 

By using the method of variation of the parameters, one can find the following solution of 
equation ([77|l in the separate variable form (f29|) 



w{r, 6) = rU Ci cos 29 + ^2 sin 29 + ^+ (78) 

C4 



^ 2A2 



(A^ sin^ 9 + cos 20)^ + \/l - A2 (sin 26* log ^2(1 - A^sin^ 6*) - 6^ cos 26* J 



The boundary conditions (fT5]) and (fT6|) necessarily imply that Ci = C3 = C4 = 0, so that 
the corresponding solution reduces to 

w{r,9) = C2r'^sm29. (79) 

where C2 depends by the amplitude of the loading (fT7|) and also by the constant associated to 
the term s = 4, that is neglected in our analysis. 

Appendix A 

In this Appendix, we demonstrate that the dynamic J-integral expression ()23p for a Mode III 
steady-state propagating crack in couple stress elastic materials is path independent. 

Considering a closed oriented path formed by two crack tip contours Fi and F2 and by the 
segments of the crack faces of length d that connect the ends of these contours, the energy 
flux integral corresponding to this entire closed path Ttot for a Mode III steady state crack in 
couple-stress materials is given by 

F{Ttot) = F(T2) - F{Ti) = V ([ [{W + T)n^-t'^n-e,w^^-fJn-(f^^]ds, (80) 

then from the definition ()23p we derive 

J{Ttot) = J{T2) - J(Fi) = (f [{W + T)n^ - t^n • e,w^^ - fi^n ■ ^J ds, (81) 

where the notation i7(Fi) and J{T2) denotes that the dynamic J-integral ()23p is evaluated 
respect to the crack tip contours Fi and F2, respectively. Applying the divergence theorem to 
the ()8ip we obtain and remembering that n^; = n • e^,, we obtain 



JF2) - J(Fi) = / V-[{W + r)e, - t^e.w,, - fi^ifj 

JAtot 



dA, (82) 
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where Atot is the area within the closed area. For the antiplane steady-state problem the strain 
elastic energy density and the kinetic energy density are given by 

W = -{aisW^x + Cr23W^y + /illV3l,a; + fJ.12V'2,x + IJ'21^1,y + IJ'22^2,y), (83) 

T=y(p^% + M,. + M,.), (84) 

the first term of the integral ()82p is the given by 

V • [{W + r)e^.] = {W + T)^^ = V^ {pW^xxW,x + J^l,xx^l,x + J^2,xx^2,x) + 

+ 2 {^l3,xW,x + CrisW^j;^; + a23,xW,y + Cr23W3,yx) + 

1 , , . 

+ 2 \l^n,x^l,x + fJ'll'Pl,xx + fJ'12,x'P2,x + fJ'12^2,xx+ (85) 

+ fJ'21,x^l,y + fJ-21^1,yx + fJ-22,x^2,y + fJ-22y^2,yx) ■ 

Taking into account the dynamic equilibrium conditions ([5]), the second term can be written as 
follows 

V • (t^e^w,^) = (V • t^) • e,w,^ + t^ ■ Vw^^ = 

= pUW^x + (o-13 + ns) W^xx + (-723 + ^23) W,yx = (86) 

= Pv'^W^ccxW^x + (O-IS + T13) W^xx + (o-23 + T23) W^yx, 

while the third 

V • ifi^ip^x) = (V • /x'^) • (^,,. + /i^ • Vip^x = 

= {Jipi - 2T23) ^l,x + {J'P2 + 2ri3) ip2,x + fJ-U^l^x + Ml2'/?2,a; + P2Wl,y + l^22^2,y = (87) 

= [Jv^i^l^xx - 2r23) Lpi^x + {Jv'^(p2,xx + 2ri3) (p2,x + l-l'Ufl,x + ^12V'2,x + fJ'2lVl,y + P22V2,y 

Substituting ()85p . ()87p and ()87p into the integral ()82p and writing ipi and 932 as functions of the 
displacement by means of relations ([I])(3^4), we obtain 

Ji'^2) - J(^l) = -^i (^13,xW,x + (^23,xW^y + - {pil,xW,yx - Pl2,xW,xx + P21,xW^yy - P-22,xW^xy) j - 



i(..........^......-............-.....) 



dA, (88) 



finally, introducing into the ([88|) expressions ([8]) and ([9]), defining the stress and couple-stress 
tensors in function of the derivatives of the displacement, we get: 

Ji^2) - J(ri) = -^r- / [{W^XXW^XXX + W^yyW^yyx " T] {w^yyxW^xX + W^xXxW^yy)) " 

^ JAtot 

- {w^xxW,xxx + W,yyW,yyx - V {w ,yyxW ,xx + W ,xxxW ,yy))] dA = 0. (89) 

We have demonstrated that for a Mode III steady state crack propagation in couple-stress 
elastic materials the difference between the energy release rate calculated considering two dif- 
ferent paths around the crack tip is zero, as a consequence we can conclude that the J-integral 
expression ([^5|) is path independent. 
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Appendix B 



In this Appendix we demonstrate that the two alternative expressions for the J-integral (j56p . 
one of which is in function of the tractions and the other of the reduced tractions, are equivalent. 
The equivalence of the two forms is demonstrated for a general three-dimensional steady state 
crack problem, and the expressions (|56|) . valid for the Mode III, are derived as a particular case. 
On the basis of general expression for the energy flux ()19p , for a steady state three-dimensional 
crack propagating in couple stress elastic materials the J-integral is defined as follows: 



J= I [(M/ + TK - t^n • u ^ + /x'^n • yj ,,] ds. 



(90) 



An alternative expression for this integral is given in function of the re duced tractions p and 
of the couple stress tractions q (JGeorgiadisl . 120031 : iGourgiotis et al.l . I2OIII ) : 



J= I [{W + T)n, 



pu^-q- cpj ds. 



(91) 



We now demonstrate that expressions ()90|) and ()9ip are equivalent. Remembering the definition 
of the reduced tractions and of the couple stress tractions ^, the following relations are derived: 



p.u, = t^n.u, + -(V/.„„xn).u. 



.T^ 



Making the difference between expressions ()90p and (|9ip. we obtain: 



(92) 



- / [{Vfinn X n) • u,x'] ds- (/i„„n • ip^,j.)ds = 

If ^ f 

^ Jt ^ Jt 



(93) 

where Sijk and £kji are elements of the Levi-Civita tensor, and the components of the rotations 
vector if have been expressed ii i functi o n of t he displacements according to kinematic compat- 
ibility conditions introduced bv iKoiterl ( 19641 ) and reported in Section [2] for the antiplane case. 
Taking into account the permutation properties of the Levi-Civita tensor elements Ekji = £jik 
and Sijk = —Ekjii the PSJ) becomes: 



If 1 /" 

If ^ f 

If If 

"n^kji I \l-^nn,i'n'j'^k,x ~r ^•n'nJ^j'^k,xi) ^^ — 77 / ^j^kji \f^nnUk,x) j ^^- V^^j 

^ Jr '^ Jv 
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Defining the vector a = ^ri.n^,x of components at = fJ-nnUk,x, the (fM|) can be rewritten as: 
- 2 / '^j^kjiak,ids = "2 / njEjikOk^ids = "2 / rijhjds, (95) 

where 6j = £jikCik,i are elements of the vector b = rota. Applying the divergence theorem, it 
follows that: 

-]- [ h-nds = -l [ {V -h) dA = -l- / [V • (V X a)] dA = 0. (96) 

2 Jr ^ Ja ^ Ja 

We have demonstrated that the diff'erence between the two alternative forms of the J-integral 
(j90p and ()9ip is zero, then expressions (j90p and (j9ip are equivalent. For a Mode III crack, using 
the definition of the out-of-plane displacement (J12p . the (|9Up becomes: 

J= f [(W + T)n^ - t^n • e,w,^ - fi^n ■ ^^^] ds. (97) 

In the same antiplane case, the following relations are derived: 

P • ^,x = P3 ■ w,x = (p • fiz)'^w ■ e^, (98) 

q • V.x = QiVi,x + q2^2,x = [(Vv')^q] • e^;, (99) 

then, substituting these expressions in equation (|9ip . we finally obtain the second form for the 
J-integral (piP : 

J = / [{W + T)n-p- e^Vw - (V(^)^q] • e^ds (100) 
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